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ABSTRACT 


Three models of an aircraft repair facility are developed 


for use in computing manpower savings achievable by scaling- 


up aircraft maintenance shops. One model is a bivariate 


Markov process model requiring a Gauss-Seidel iterative 


algorithm for its solution. The other two 
Repairman problem and the M/M/S queue, the 
can be found in most introductory texts on 


The simple models are compared to the more 


are the simple 
solutions to which 
Stochastic models. 


refined bivariate 


model in regard to predictive accuracy. The simple models 


are used to compute the manpower savings achievable in one 


scaling-up scheme, and predict a 20% savings in maintenance 


personnel. 
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TABLE OF SYMBOLS 


K total number of shops at the repair facility. 
M total number of aircraft at the repair facility. 
B as a subscript denotes back-shop activity. Also 


used as a random variable indicating the number 
of aircraft down at a shop for back-shop repair. 


F as a subscript denotes flight-line activity. 
Also used as a random variable indicating the 
number of aircraft down at a shop for flight-line 


repair. 

r arrival or failure rate 

R the set which characterizes the repair shops 
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portion denoted by the subscript s 

u repair or service rate 

P the set of performance measures which characterize 
the repair facility 

E(arg) the expected value of the argument 


Var(arg) the variance of the argument 
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indicated by the subscript i 

WwW. waiting time at the shop indicated by the 
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Fa 


‘Ba 


the number that the random variable B assumes 
the kth iteration value of P ey 


the long-run probability of f aircraft down for 
flight-line repair 


the long-run probability of b aircraft down for 
back-shop repair 


the long-run probability of n aircraft down for 
repair (n = £ + b) 


the expected number of aircraft in the system 
the average arrival rate at the flight-line 
the average arrival rate at the back-shop 


the average arrival rate into the shop, 
flight-line and back-shop both 
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I. INTRODUCTION 


It is quite a three-pipe problem. 
-Sir Arthur Conan Doyle, The Red-headed League 

A. BACKGROUND 

Under the auspices of the United States Air Force Project 
Rand, the RAND Corporation is currently studying ways to 
improve the aircraft maintenance posture of the United States 
Air Force. In the past, RAND's Logistics Studies Program 
has been primarily concerned with the improvement of particu- 
lar logistics functions. They have since departed from this 
basically micro approach to logistics problems for several 
reasons. First, the Air Force has over the years developed 
an in-house capability to analyze micro-oriented logistics 
problems, so that RAND assistance in this area is no longer 
required to the extent that it was in the past. Second, 
with the Air Force eeuietehiee analysts now concentrating on 
the micro-oriented problems, RAND's analysts are free to 
concentrate their efforts on macro-oriented approaches to 
the total logistics system. "If the Air Force can take care 
of the trees, then let us study the forest," is a good, 
although casual, expression of the present RAND approach 
to logistics problems. This ivory tower attitude may result 
in the introduction of new logistics structures in the future 


that will dramatically reduce support costs. 


12 





This shift in emphasis from a micro approach to logis- 
tics function, to a concentration on macro-oriented problems, 
is a fairly recent development. In July, 1973 a project 
called ANALOGS 1980s (Analysis of Logistics Structures for 
the 1980s) was formally initiated as a part of RAND's 
Logistics Studies Program. Sponsored by Lieutenant General 
William W. Snavely, Deputy Chief of Staff, Systems and 
Logistics, Headquarters, USAF, the project has as its explicit 
goal the exploration of logistics structures that could 
dramatically reduce support costs without degrading mission 
capability. Later, on September 9, 1974, the Deputy Chief 
of Staff, Systems and Logistics, provided further impetus 
to the macro approach when he issued a Program Management 
Directive authorizing a Maintenance Posture Improvement 
Program (MPIP) with two broad objectives: 

1. To reduce costs of maintenance in peacetime and 
Wartime contingencies. 
2. To increase the effectiveness of maintenance in 
Support of operational missions. 
Information on ANALOGS 80s and on work to support the MPIP 
has been provided by S. Drezner and I. Cohen of RAND, to 
whom the author is grateful. 

As a result of RAND work done in support of MPIP and 
the ANALOGS 80s project, there is a current interest in 
two organizational structures which attempt to reduce support 


costs by scaling-up maintenance activities, both of which 
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represent a departure from the organizational structures 
presently employed by the USAF. The older structure of the 
two, called Queen Bee, specifies the centralization of cer- 
tain intermediate level aircraft maintenance at large central 
bases called Queen Bees. The second structure, called the 
Reallocation of Activities Alternative (RAA), grew out of 
the ANALOGS 80s project analysis. RAA specifies the decen- 
tralization of combat-mission-oriented activities at Combat 
Mission Bases (CMBs), and the centralization of logistics 
and training activities at Support Mission Bases (SMBs). 
Costs and payoffs calculations for each of the two struc- 
tures are needed before a decision to implement is made 
regarding either structure. In late September, 1974, work 
began in earnest on the development of models which could 


be used to assess the effectiveness of such configuraticns. 


B. MOTIVATION 

At the outset of the modeling effort, it was felt that 
the models developed could be used in a broader context than 
the RAND study. The Air Force is not the only branch of the 
armed forces with an aircraft logistics problem; the Navy 
has one too. Additionally, there were similarities between 
the structures employed by the two branches, so that the 
models developed for the Air Force could be applied to Navy 
problems with little or no modification. There was, thus, 
strong motivation for a thesis focusing on the modeling of 


aircraft logistics structures. 
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C. USES 


imeiclySls Of euLrent Structure 


Once models for the current structure have been 
developed, they could provide answers to several cogent 
questions. For example: Is the present manpower allocation 
within the current structure reasonably satisfactory, if 
not optimal. How might the manpower be reallocated to 
improve Peeaaneser What degradation in readiness can be 
expected for various reductions in manpower? If manpower 
reductions are necessary, where will the effects of the 
cuts degrade readiness the least? 

Most importantly, models for the current structure 
also provide a basis for the comparison of the current 
Structure with a proposed alternative structure. They can 
also be used to analyze the effects of increasing the number 
of aircraft that the current maintenance structure must 
Support. Alternatively, they can be used to indicate 
acceptable manpower cuts in the event that the number of 
aircraft and/or the amount of flying is decreased. Provided 
that reliability estimates of a new proposed weapon system 
were available, the mcedels would be useful in predicting 
manpower costs in support of a new weapons system. It may 
be that an alternative support system may be better suited 
for the proposed weapon system than the current support 
system. The models would help resolve problems such as these. 

2. Analysis of Alternative Structures 
As stated previously, before decisions are made 


regarding alternative structures, the current system should 


ibs 





be modeled to provide a basis for comparison of the current 
structure with the proposed structure. In the RAND study, 

the current structure was modeled first. This provided 
measures of performance for the system. Then the alternative 
structure was modeled. The model of the alternative struc- 
ture was used to show that measures of performance on par 
with those attainable with the current structure were achieva- 
ble with the proposed structure with fewer personnel. 

There is a trap to watch for in this kind of com- 
parative analysis. If the manning of the current support 
system is not optimal, then the comparison with the alterna- 
tive structure, assumed operating optimally, injects a bias 
in favor of the alternative structure into the analysis. 

It may be that the reallocation of manpower within the current 
structure will yield measures of performance that the 
alternative structure cannot satisfy with less personnel. 

3. Application to Other Support Structures 

It was stated earlier that there were similarities 
between Air Force and Navy aircraft support structures. 

There are also similarities between aircraft support struc- 
tures and ship support structures. If there are economies 

attainable by scaling-up certain logistics functions within 
the aviation community, then similar economies may exist in 


the surface and submarine communities. 


D. SCOPE 
A logistics support structure has not been defined up 


to this point. A logistics support structure is defined as 
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the organization, training, quantification, and geographical 
positioning, where each 1s applicable, of military and 
civilian personnel, tools, equipment, armament, spare parts, 
fuel and physical plant facilities with the explicit pur- 
pose of maintaining a weapon system in a satisfactory state 
of readiness to perform peacetime missions and meet wartime 
or combat contingencies. This nebulous definition is pro- 
vided to point out the magnitude of the task confronting 

the analyst who aspires to model the structure in toto. 

It was decided to look at a small but significant part 
of the total structure, and see if there were economies to 
be derived from the scaling-up scheme proposed in the Queen 
Bee and RAA structures. To restrict the analysis to the 
organization of repair personnel and calculate manpower 
savings achievable with the alternate structures appeared 
to be the best approach. This approach was appealing for 
Several reasons. First, support personnel costs are a 
Swemericant contributor to the total lifetime cost of a 
weapon system; a significant reducticn here can dramatically 
reduce the cost of a weapon system. Second, repair organi- 
zations can be modelled as queueing systems. There are 
standard calculations of measures of congestion associated 
with queues which makes side-by-side comparison of two 
different repair organizations relatively simple. Third, 
many other costs associated with the logistics structure Age 
pend in large part on the numbers of repair personnel and on 


their organization and positioning within the logistics structure. 
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II. DESCRIPTION OF CURRENT REPAIR ORGANIZATION 


Why is this thus? What is the reason of this thusness? 
~Artemus Ward, Artemus Ward's Lecture 

The current aircraft repair organization, hereafter 
referred to as an aircraft repair facility, consists of K 
shops, designated Shop 1, Shop 2, and so on, up to Shop K. 
There are two basic categories of maintenance tasks performed 
at each shop: flight-line repairs and back-shop repairs. 
The flight-line repair category consists of all those repairs 
that can be made on board the aircraft at the flight-line. 
The back-shop repair category consists of all other types 
of repairs. The repair facility supports an operational 
unit which operates a fixed number of aircraft, equal to M. 

Figure 1 is a schematic diagram of the operational unit 
and a repair facility composed of six shops which shows how 
the facility and the operational unit interface. When a 
mission is scheduled, the required number of aircraft are 
removed from the pool and preflight maintenance is conducted. 
When preflight maintenance is satisfactorily completed, the 
aircraft conduct the missicn, return to base where post- 
flight maintenance is conducted. Once this is completed, 
they are returned to the pool. During stand-down periods, 
other maintenance is conducted in addition to the preflight 
and postflight maintenance. This maintenance is referred 


to as scheduled maintenance. The portion of the diagrammed 
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repair shop with enclosed capital B represents the back-shop 
activities of the repair shop, and the portion with circled 
numeral and capital F indicates the flight-line activities. 
Small f and small b denote flight-line and back-shop failures. 
When a failure occurs, the cognizant shop dispatches a 
repairman to the aircraft where the failure occurred. If 

it can be repaired onboard the aircraft it 1s designated a 
flight-line failure. If not, then it is considered a back- 
shop failure. 

Back=-shop or flight-line casualties can occur during 
missions, preflight, postflight or scheduled preventive 
maintenance. ‘We assume that casualties do not occur to 
aircraft that are in the pool or already down for repair. 

It mgs also assumed that the tempo of flight operations is 
such that the number of aircraft in the pool is very small 
compared to the number of aircraft available to the opera- 
tional command, so that the total number of aircraft which 
can fail at any given time is approximately equal to the 
total number available, M, minus the number that are down 
for repair. 

A repair facility consisting of K shops, each capable 
of performing two categories of maintenance tasks implies 
that there are 2K possible failures that the facility can 
handle. It is assumed that the 2K possible failures exhaust 
all of the types of failures that can occur on the aircraft. 


It is further assumed that all aircraft at the operational 





command exhibit identical failure characteristics; all 
aircraft fail into the flight-line portion of Shop 1 at one 
constant rate, and into the flight-line of Shop 2 ata 
different constant rate, and so on, for all K shops. Back- 
shop failures exhibit the same property. 

Each shop at the facility is characterized by the number 
of repairmen in the shop, and two repair rates, one for the 
flight-line repairs and the other for back-shop repairs. The 
repairmen are assumed to be equally capable at any task 
confronting the shop. This makes it possible to characterize 
the shop by specifying the number of repairmen, the flight- 
line repair rate, and the back-shop repair rate. Repairmen 
at one shop cannot assist repairmen at a different shop, so 
that if failures occur when all repairmen are busy at one 
shop a queue begins to form. 

A priority system within each shop is in effect to insure 
that repairmen are assigned to flight-line repairs before 
commencing work on back-shop repairs. Furthermore, every 
repair task can be accomplished by one and only one repairman. 

Let F denote the set of 2K different failure rates which 
Characterize each of the M aircraft; then 
r r 


r hee Jee CZ.) 


= { FK’ “BK 


APL! pales r2'“p2! ers 


where the alphanumeric subscript indicates the type of 


failure, F and B denoting flight-line and back-shop 


Zr 





respectively, and the numeral denoting the cognizant repair 
shop number. 


Let R denote the set which characterizes the repair shops. 


Then, 
lo  \ 
S, s.\ S S, oe 
R = Yny g Heo g a aie g tle | a t Pax CZ 2) 
Cr lsd Sify) 3.8) PBK/ | 


where Sy is the number of repairmen, is the flight-line 


“Toa 
repair rate, and WRG is the back-shop repair rate at Shop 1. 
Another set of interest contains measures of congestion 


at the shops in the facility; it is designated P and is 


defined as follows: 


oS {P,, Pays Pai ee ley Zep Sighe ety (2753) 
where 
E(N, ) E(Nas) E(N,,) 
Var (N, ) Var (N,,;) Var (N,;) 
q E(Q.) 7 E(Q.,.) E(Q..) 
2 = , Pas = Ba , and ae = ee 
Var (Q;) Var (Q,,;) Var (Q,;) 
E(W;) E(W.,;) E(W,;) 
E(D;) | E(D,,, ) | E(D,;) 
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The subscript i refers to the ith shop, and the subscripts 
Fi and Bi refer to the flight-line and back-shop activities 
of the ith shop. The variable N is the number in the system, 
and the variable Q is the length of the queue. W is the 
Waiting time, and Dis the delay time. Waiting time is 
defined as the total time the aircraft is down for repair. 
Delay time is the time spent waiting for repairs to begin. 
Figure 2 illustrates the interaction of the elements of 
set F with the set R, resulting in the performance measures 
in set P. Note the feedback loops out of R; these represent 
the interaction of the shops with each other. This inter- 
action is not well understood at present, and models of the 
present facility developed on the follcwing pages ignore 
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III. A BIVARIATE MODEL 


How happy could I be with either, 
Were t'other dear charmer away! 
But while ye thus tease me together, 
To neither a word will I Say. 
-~John Gay, The Beggar's Opera 

A bivariate model is developed in order to obtain values 
for the measures of congestion contained in the set P. The 
basic approach in modeling the system with a bivariate 
model is to treat each shop independently, and model the 
facility as K different queueing systems, each queueing 
system maintaining M aircraft. Figure 3 illustrates the 
approach. 

Consider the queueing system defined by the interaction 
of Shop N with M aircraft. Assume that repair times and 
interarrival times are exponentially distributed. The state 
of the system at any time t can be represented by the ordered 
pair (F(t) ,B(t)), where F(t) is the number of aircraft down 
for flight-line repairs, and B(t) is the number of aircraft 
down for back-shop repairs at Shop N. If flight operations 
and maintenance are conducted throughout the 24 hour day, 
so that the process is not dependent on time, the bivariate 
model belongs to a class of probability models known as 
continuous-time Markov chains. These probability models 
are characterized by the fact that 

1. They can be completely described at any time by 


specifying their state at that time, and 
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2. the time until a transition is made from one state 
to a neighboring state is an exponentially distributed 
random variable. 


Let 


P ep (t) = P{F(t) = f, B(t) =b]). (35.1) 


Concentrating on long-run probabilities, of more interest 


is the quantity 


P eb = lim P oy (t) (332) 


Hee te (t) = £,.. B(t) 


t->0o 


b). oo) 


Pep is the long-run or steady-state probability of exactly 
£ aircraft down for repair at the flight-line, and b aircraft 
down for repair at the back-shop of the shop under considera- 
Cwen . P ey is also the long-run proportion of time that there 
are exactly £ and b aircraft down for repairs at the flight- 
line and back-shop respectively. 

For a queueing system consisting of one shop servicing 
M aircraft, the possible states of the system are (f,b), 
Pee, 2;2,2--, M, b= 1,2,3,..., M, such that £ +b < M. 
The probability distribution of interest is the steady-state 
probabilities associated with the possible states of the 


system 
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The number of possible states is a quadratic function of the 
number of aircraft serviced by the shop. Denoting the 


number of possible states by Noo the relationship is 


OM ee eee eee) 
ir 5 : (3.4) 
The state space for a system with a large number of aircraft 


can be guite large, as Table I illustrates. 


LABLE & 


Number of Bivariate States 
Vs. Number of Aircraft 


M N. 

3 

Zu 

10 66 

25 Soe 

50 vs26 

f° Z926 

100 Silieyt 
200 20301 


The general principle which enables one to determine 


the state probabilities is derived from a rate-equality 
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principle which applies to continuous-time Markov chains. 
The principle applied to the bivariate process states that 
for each possible state (f,b), the rate at which the process 
enters the state (f,b) equals the rate at which it leaves 
the state (f,b). The application of this principle to each 
of the No possible states of the system will yield a set 

of balance equations (N. in all) which must be solved for 
the set of state probabilities {Pept 

Figure 4 is a transition diagram for a hypothetical 
shop consisting of three repairmen with repair responsi- 
bilities for five aircraft. This shop represents the repair 
facility described in Chapter II. The circles represent 
the possible states of the system. Arrows indicate the 
possible transitions between states, and the labels on the 
arrows represent the rates at which the process moves from 
one state to another. The transition diagram simplifies the 
task of setting up the balance equations. 

Balance Equations for a system of M aircraft failing 
into a shop with S repairmen are presented in Figure 5. 
These can be verified by referring to the transition diagram 
of Figure 4. 

The solution of the system of balance equations can 
proceed algebraically utilizing a successive substitution 
scheme for small systems, but for large systems the solution 
becomes intractable. For this reason, an iteration algo- 


rithm was invoked so that {Pept for systems of practical 
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size could be found. The algorithm used was the Gauss- 
Seidel iterative method described by Varga in [Ref. 1] and 
by Young and Gregory in [Ref. 2]. The application of itera- 
tive algorithms to Markov-type problems is an area of active 
research at the Naval Postgraduate School, as it is particu- 
larly well suited for the solution of large systems of 
balance equations. 

The algorithm as it applies to the problem at hand is 
set forth here without mathematical underpinnings. First, 
all balance equations are rearranged so that P en 1s expressed 
as a function of eS P £+1)b! Pe (b-1) ’ and PE (D+1)° 
This 1s easily done, as reference to equations (3.5) and 
(3.6) will show. Second, assume values for all Pep: One 
scheme would be to assign P ey = 1/N,, EOmmadin & and» 5b in 
the state space. Here, N. is the number of possible states. 
The closer the initial assignment of values of P ey is to the 
solution, the faster the algorithm wili converge. Third, 
denote the assumed values of Peper Bea for the zeroth 
iteration. Fourth, establish a scheme for incrementing sub- 
Scripts £f and b. The scheme used in the problem under con- 
Sideration was to start at P and increment £f values by one 


00 


until P was reached, then b was incremented by one, and 


MO 
f values were then incremented by one, starting at f = 0 
again. Thus, after Po was reached, the incrementing scheme 
called for Pol! Poa Poi en a EM-D1 in that order. The 


subscript b is again incremented by one and f values 
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incremented by one, so that £f = 0,1,2,..., M-2. This process 
is continued until all Pep in the state space have been 
reached. This constitutes one iteration. So one iteration 
will consist of N. increments. Unless we do something at 
each increment, of course, the algorithm would be rather 
pointless, so on to the remaining steps of the algorithm. 
Denote the value of Pep computed at the nth iteration 
_—". Recall that all balance equations were rearranged 
so that P ey equals some function of 2 BTN Se P (£41)b! 
Pe (p-1)’ Pe (D+1)° Bach balance equation will generate a 
different function, but to keep the notation simple, denote 
that function generated by the balance equation under con- 


Sideration g, where it is understood that the gs will be 


different for each equation. We have 


Pep = TP ce-aypt Pcetiyy’ Pe(b-1)" Pe(pt1)?* cel) 
Define 
(n) _ (n-1) (na) Goi IL) Coll) 
ms 8 (eam) | f+1)b 2S jolly C2 se dere 
(3.8) 


as the nth iteration approximation of Pep: Note that it 
is computed from the (n-1)th iteration approximations of 


and P 


Mee t)b’ § (£41)b’ *£(b-1)" £(b+1)° 
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The nth iteration values of Pep are computed as follows: 
ime For the first P eb in the incrementing scheme, compute 


P (n) in accordance with equation (3.8). 


fb 
2. Immediately change Berea. aa) Dene For example, if 


P (1) is being computed, then immediately after the step 


00 
where a value for ee a is obtained, change Bea’ to that 


value. Now, when the incrementing scheme calls for the 


a) (0) 


computation of a P eb involving a Poo , the algorithm will 


use the most recent approximation of Poo? which is the 


(1) 
00 ° 


3. Continue through the incrementing scheme, doing steps 


value computed for P 


1 and 2, until the algorithm has incremented through all 
possible states. 
4. Make the following computation: 
M M-b 


-_- =) Bee (3.9) 
b=0 £=0 


— by dividing each one by 
(n) , 


the value of S obtained in step 4 above. The P ey s now 


5. Normalize all of the P 


Sum to one. 


(n) 
fb 


the computation of ye during the (n+l1)th iteration, 


fee ane values of the P just obtained are used in 


which proceeds exactly the same as the nth, following steps 


imetnrough 5. 
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Eventually, the algorithm will converge on the solution 


i The bivariate model 


of the system of balance equations. 
balance equations were solved in the above manner. The 
average number of iterations required for convergence was 
approximately forty-to-fifty, with one system converging 
in seventeen iterations. The end result of the algorithm 
is a solution for {Ppp} which was used to calculate the 
measures of congestion contained in the set P (equation 
3) ) . 


The calculation of the elements of P is straight forward. 


Each component of the vector elements P., Bene 


and Pai is 

a standard queueing system calculation. The computation 

of each component proceeds without difficulty, once the set 
{Pp} is known. 

Define {Pe} as the univariate probability distribution 
of the flight-line portion of the queueing system under 


consideration. Then 
{Pe} ieee ieee t = O,1,.2,..65M , (3% 10) 


where 


the fact that the Gauss-Seidel iterative method 
described converges has been established by J.P. Lehoczky 
(private communication). 
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Mar 


P. = P{[F =f] = & P (3.41) 
£ b=0 fb 
Similarly, for the back~-shop, 
Eee Pies, b> = O0,1,2;,...,M ; C3282) 
where 
M-b 
ee P([B = bj] = ae Pep - (372153) 


The univariate probability distribution of the number 
of aircraft in the queueing system in the aggregate 


(flight-line and back-shop) is defined as Pat where 


Po ie Wines 0 = 0,1,2,6..,;M-, (3.14) 
with 
n 
P = P[N = n] z ee. (ogene (3245) 
£=0 
or 
n 
12) Sta 2, a Saree : C3216) 
n b=0 (n-b) ,b 
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The expected number of aircraft down for flight-line 
repair at the queueing system under consideration (shop 


1) is 
M 
E(N,, | = a fe Pe (3e17) 


where Pr, is the value obtained from equation (3.11). 


Similarly for the back-shop, 


| 4 = 


E(N, 1 = 


b P (Om Le) 
b b 


if 


where Py is the value obtained from equation (3.13). The 
expected number of aircraft down in the system as a whole 
is 

E(N,] = n P (2 JUS) 
where aa is obtained from either equation (3.15) or 
equation (3.16). 

The expected number of aircraft in the flight-line queue 
is 
M 
Bl@=i = ££ (£-—- S.) P : (S20) 
Fi f=S.4] aL fe 
a 

The expected number of aircraft in the back-shop queue is 


obtained after some head scratching. When the number of 
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aircraft that are down on the flight-line is less than or 
equal to the number of repairmen at the shop (f < S35), 
f repairmen will be working at the flight-line, with the 
remaining S; - £ men available for back-shop repairs. As 
long as the number of aircraft down for back-shop repairs, 
b, is less than or equal to Ss; f, no back-shop queue will 
exist, Since there is an adequate number of repairmen 
available to handle back-shop tasks. This implies that the 
@meue at the back-shop exists when b > S. = 3Ee a dan eet 
queue will contain (f + b - S;) aircraft. Thus the expected 
number in queue at the back-shop when 0 < f S84 is 

Si M+ 


oy o (b+f£-S, ) 18 


£=0 b=S,-£+1 fb . 


When f 2S. all repair personnel are tied-up with flight- 
line repairs, so that every aircraft down for back-shop 
repair is in queue, awaiting service. This yields 

M-1 M-f 


2 toy 12 


f=S$.+1 b=1 Eb 


Adding the two expressions together gives the expected 
number in the back-shop queue: 
Si M-£ 


BO Fd y (bo + £ - S.) P 
a £=0 b=S,+1-f i ~ fb 


Mek yea 
(3220) 





The expected number of aircraft in queue for the system 
as a whole is Similar to the back-shop queue in regard to 
derivation. When f <S;, the only aircraft in queue are 
those in the back-shop queue, Since there are enough 
repairmen to handle flight-line repairs. Now, when f > Ss 
queues can exist at both the flight-line and the back-shop. 


The number of aircraft in queue will be (f +b - S;). This 


yields 
Si M-£ 
E(Q.J] = £ r (b + £ - S.) P 
+ £50 b=S, +1-£ i” “fb 
M M-£ 
ny Min f£ Ses P. (3920) 
f=S +1 b=0 iv fb 


The variances of the number of aircraft in the flight- 
line, back-shop, and total system are VariINnsl, VarIN,.], 


and Var[N,], where 


: 2 2 
Var IN, ] = EIN, ] (EIN,; 1) (3.23) 
M , M 2 
oooh. - Prt r,) , (3.24) 
aa ash 
Var{N..] = EIN..*] - (E(N..])° (3.25) 
Bil BL BL _ 
M M 2 
= 7 ep 2 =) 8 eee (3.26) 
eal b=1 
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and 


_ 2 2 
Var [N, ] = E(N, ] (E(N,]) a2 7) 
M 5 M 2 
> on P =e A Ro oe . (3.28) 
n=l = n=1 ss 


The variance of the number of aircraft in queue at the 


flight-line is 


2 
2 
varlQ.,1 = EQ, 1 (On 1) (3.29) 
M . M 2 
= f (ee CiemE. Felis GE-sa Pe). 
f=S +1 = f f=S +1 it 


(3330) 


The variance of the number of aircraft in queue at the 


back-shop is 


2 
2 2 _ 
Var lQ,;] = E[Q,; ] (E[Q,.]) (32a) 
So; Met : 
= 2 ye (b+f£-S. ) Pep 
£=0 b=S.+1-f 
M=-l M-f£ 5 
+ YF RD ©: P ep 
f=S.+1 b=1 
1 
2 
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The variance of the number of aircraft in queue at 


the system as a whole is 


2 


_— 2 = 
var[Q,] = E[Q,°] - (E[Q,}) (3.33) 
Si M-£ 2 
=p Gael) lags 
£=0 b=S,+1-f 1 Eb 
M M-£ 2 
+ fF ys (f+b-S.) P 
f=S.+1 b=0 eS 
1 
2 
- (E{Q,1) . (3.34) 


The calculation of the expected waiting times E(Woid, 


BE [W i! and E[W.] is based on the queueing relationship 


B 
"IL = AW", where Lis the expected number of customers in 

the system, » is the average arrival rate, and W is the 
expected waiting time. Denoting the average arrival rates 

at the flight-line and back-shop by Pa and Ba! respectively, 
and the average arrival rate to the total system by hat 


the queueing relationship is restated in terms of notation 


applicable to the problem at hand. Thus, we have 


E (Np, | = pa E(w], (3.35) 


E(N,.1] ee E(wB.], (3.36) 
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EIN,] = A, ElW,]. (3.37) 


Given that the state of the system is (f,b) there are 
exactly (M - f - b) aircraft that can fail into the system. 
They each have individual failure rates hei? so that the 
failure rate into the system in state (f,b) is (M- £ - D) Ag. - 
The proportion of time that the system is in state (f,b) 
is Pep: Summing up the failure rate into the system at 
state (f,b) times the proportion of the time that the system 


is in state (f,b) over all possible states yields the average 


failure rate for the system under consideration. Thus, 


M-f OM 
r = > & (M-f-b) A_. P (3.37) 
Fa eon e=0 Pi.  £D 
M-£ M M-f M 
= hk. ££ EMP ee se (Ce + b) 
Fi b=0 f=0 fb Fi b=0 £=0 fb 
M-f OM M-f OM M-f OM 
=i TUR Ss rE fP. = F Eb PP.) 
Pi p=0 f=0 = peo t=0 ©? peo t=9 © 
= hos [M - E(N,; 1 - E(N,,1). (3.38) 
Ey. Z 3.39 
ee Ap; %[M - EIN,]] . ( 
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Similarly, 


Apa = a [M - E(N, J] ° (3.40) 


Derivation of the overall system average arrival rate 


yields 
M-f£f M 
N- = con Pls as = i = b) Ons ar hpi? P en (3.41) 
M-£ M 
n= eee ae ri * tp 
M-£ M 
= Pa “ Ba (342) 


which is no surprise. 
The calculation of delay times is based on the queueing 
relationship Q = XD. Stated in terms of applicable notation 


the delay times are: 


E[Q,,] r E[D,.] ’ (3243) 


Fa 


r 


E(Q,;] E[D,,] ; (3.44) 


Ba 





and 
E[Q,] = rs E[D; ] : C3..45) 


The bivariate model is an approximate model of the 
current aircraft repair facility described in Chapter II. 
The purpose of the development of the model is to provide 
a basis of comparison for scaled-up versions such as the RAA 
and the Queen Bee structures. Two other approximate models, 
called decomposition models, were developed to assist in 
modelling the scaled-up structures. A discussion of these 


models is presented in the next chapter. 
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IV. TWO DECOMPOSITION MODELS 


He was in logic a great critic, 

Profoundly skilled in analytic. 

He could distinguish, and divide ... 

-Samuel Butler, Hudibras 

Initially conceived to provide a rough-and-ready model 
of the current repair facilities, the decomposition models 
were soon abandoned in favor of the more refined bivariate 
model described in the last chapter. However, if scaling-up 
is to occur, the current repair shops must be physically 
decomposed in some manner, aS repair personnel will be taken 
from shops as currently organized and reassigned into Queen 
Bee or RAA Support Mission Bases. One scheme is to maintain 
sufficient personnel at each shop in the facility to handle 
flight-line repairs only, and reassign the remainder to an 
aggregated facility with exclusive responsibility for back- 
shop repairs. Men would be drawn from several bases for the 
aggregated back-shop. The idea is that the pooled back~shop 
resources from the bases would provide more than sufficient 
manpower to give adequate back-shop support. The excess 
manpower would then represent the manpower savings achieved 
by reorganization into the new structure. 

Part of the problem is to determine how many men to 
leave at the flight-line facilities. The bivariate model 
helps in this, since it provides a measure of flight-line 


performance at the shops represented by the vectors Pras 
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i = 1,2,3,..-,K. In addition, the bivariate model gives 
measures of back-shop performance in the vectors Pair 

m= 1,2,3,.--,K. A decomposition model fits in nicely at 
this juncture in the overall problem. With a model that 
yields measures of performance for various allocations of 
shop personnel between exclusive flight-line duties and 
exclusive back-shop duties, all that is needed is to find 
that allocation that meets present flight-line performance 
criteria, and the numbers available for the pooled back-shop 
are in hand. The back-shop can then be modeled, utilizing 
one of the simpler queueing models, to find the number of 
personnel required to provide comparable support to that 
provided by the current facility, and the manpower not needed 
represents the savings acquired by reorganization. 

Figure 6 best describes the decomposition approach. 
Again, independence between shops and exponentially distri- 
buted repair and interarrival times are assumed. A further 
assumption is that flight-line and back-shop activities within 
each shop are independent of each other. The S; men asSigned 


to shop i are allocated so that S men conduct flight-line 


Fi 


repairs exclusively, with the remaining S working in the 


Ba 
back-shop. 

Figure 7 is a transition diagram of the process. At 
each shop, two independent processes operate simultaneously. 


Each process is identical to that described by Gaver and 


Thompson in [Ref. 3] as the Finite Arrival Source, or 
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Repairman problem. The state of each process can be repre- 
sented by (F(t)), or (B(t)) as applicable. F(t) is the 
number down for repairs at the flight-line, and B(t) is the 


number down for back-shop repairs. Let 
P_(t) = P([F(t) = fj] . (4.1) 


Again, the long-run probabilities are the sought quantities, 


i.e. 
Py = lim [P_(t) ] 
to 
= lim P[Fft) = £] , (4.2) 
{->00 
and 
P. = lim P[B(t) = b] . C27. 3) 
£->00 


Balance equations for the queueing system represented 


by flight-line of Shop i are: 
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State Rate Process Leaves = Rate Process Enters 
0 Mp Po = Up Py (4.4) 
1 [(M-1) A. + Hp, IP, = M34P a ZU; P5 (4.5) 
2 [(M-2) A, + 2u,,; JP, = (M1)A,;P] ut 3 2s Ps (4.6) 
ee. £ < Say [(M-f) A... + fips IP, = (M-f+1) 0,;Pe_) 
+ (£41) U.P ei) (4.7). 
eet <M [(M-f)0,, + S.iup, IP, = (So) Aa 
+ Spite ey (4.8) 
: Sri tei Py [are tely .. oe! 


Replacing the term Mai Po on the left side of the balance 
equation for state 1 with its equivalent Ups Py obtained from 
the balance equation for state 0, the balance equation for 
State 1 becomes 

P 


(M-1) A (4-10) 


pit Upy Py = py Py t+ 2p Po 
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Clearing terms , 
(M-1) 25 Py = 2Ups Po (4.11) 


Replacing the expression (M-1) ABP) on the left side of the 
balance equation for state 2? with its equivalent just found 
in equation (4.11) 2upiPs , and clearing terms, the balance 


equation for state 2 becomes 


This process continued down the states from state 3 to state 


M yields the following set of balance equations: 


A-P_ - Mea Peay ; Cra lel (4.13) 
where 
Ae = (M-£) 0,5 ; 0 < £ < M-1 (4.14) 
Ue = Pups ’ O<f < Sri (4.15 
= Sri lp ; Sri S 8 <M ; (4.16) 
The solution to this set of balance equations is 
ED Ae 
2. = 2 a (eer nC (4.17) 
fe 0 j=1 ue aa 


ale 





Po is obtained from the requirement that 
M 
” Pe = J]. (4.18) 
£=0 
Thus, 
M M 
ie + t 6 6?P = Me Iz = J] 
f=1 * f=0 * 
M 2 Na 
one LP TI a = We 5 
ei jel ag 
and 
tl te Na ] 
ei l + f TI 1 | = (4.19) 
f=llLj=1 +3 | 


The solution of the balance equations obtained by 
substituting expressions for Ae and Ue obtained from equations 
(4.14), (4.15) and (4.16) into equations (4.17) and (4.19) 
yields the following solution for the flight-line balance 


equations: 


Po = I 
"Fi ag [Fea ME a M! wi \* 
1+ f | 2 
f=), weir | aut | f=S,,,+1 STS eye era || oa 
(4.20) 
Pe = GREET pee * "0 aa ae eile 


Hey 
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ve 
1 
Pe = 2. aes ise tere M 
(M-f)! S,.!S,, 0 "Fi Vey 
(Aso 


The solution for the back-shop balance equations is 
obtained employing the same method used to solve the flight- 


line balance equations. The solution is 


a ali 
Po = ’ 
So 1) 
ee, 7 M! ‘Bil : M! rei \P 
bay (PES) tb | ue b=S_.412 (Mb) is .is .P pi | Yai 
Bi B1. B21 
(4.23) 
> = is “Bi P ? 2% as (4.24) 
b (M-b) !b! Uns 0 7 ere : 
i 
and 
b 
M! Bi 
ee | BL Pa cee ab <M (4,25) 
8 (M-b) 1S, 1S.” Spi | Ua 0 oe > 
Calculation of the elements of P... and P.. were made 


Fi Bi 


utilizing equations from Chapter III where appropriate. For 
the expected number of aircraft down, equations (3.17), (3.18), 
and (3.19) were used. The expected number of aircraft in 


queue at the flight~-line was calculated from equation (3.20). 


4 
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The expected number of aircraft in the back-shop queue was 


calculated using 


M 
E(Q..) = y (b-S.) P (4.26) 
Bi b=S, +1 aL b 


The variance of the number of aircraft in the flight-line 

and back-shop were calculated using equations (3.24) and 
(3.26), respectively. The variance of the number of aircraft 
in flight-line queue was computed from equation (3.30). The 
variance of the back-shop queue was computed from 

M 


var(Q,;) = (b-S,)* P, - (EQ,,])* . (4.27) 


b=S +1 
Calculation of the expected waiting times and expected 
delay times were made utilizing equations (3.35), (3.36), 


(3.43), and (3.44) with 


en Ne E = EIN; 3] (4.28) 


and 


Nee Ans | M : EIN, 53] (4.29) 


Another decomposition model assumes a Poisson arrival 
process, with constant arrival rate Mas og Me; as applicable. 
Although not as accurate as the repairman model, it can be 


applied to systems with a large number of aircraft and 


34 





sufficiently manned so that the number of aircraft operational 
is kept at high levels. 

Again the problem has been reduced to the repeated 
solution of a simple model, this time the Poisson-Arrival, 
Parallel-server, Exponential Service Time System, or M/M/S 
Queue described by Gaver and Thompson in Ref. 4. The long- 


run probability distributions {P,} and {P,} can be found, 
£ b 














miei Zing 
p MX f 
= ff Fi 
Pe = =r |= yl © Si ESI (4.30) 
F1 
and 
P MA. fe 
= ore — oe <e (4.31) 
4 ae : 
337 Jee Fi 
where 
MA. -1 
Sry Mi f MA Fi S = 
q i Fi “By 1 FilPi 
PEO 7 » tr | Lig ar Care MA 
f=0 ~° Fi Fi Fi’ _ Fi 
Spi lps 


(4 32) 
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For the back-shop, use 











peor. pal : 
PB, = rl Crees. «= (ares) 
Bi 
b 
Pp MA. 
b0 Bl 
= So =- bp (4.34) 
b-Spi ile = / Bl — 
533! Spi Bi 
where 
S_, b SBi MBA 7 
7 Bi, /MABy MN 1 "Bi Bi 
Pho } SE u oa ar S..! M 
b=0 ~° Bl Bi Bi ea Bl 
q "pi lpi 
CARS 5) 


The expected number of aircraft down for repairs at the 


flight-line and back shop are 











S MARS 
Mi ma_,\ Fi aD a. 
E(N,;] f+ +P. (= — er sy, (4,36) 
| Mpa i) Fi [2 a Fl 
Spy bey 
and 
an MAB i 
— MARG See il "3a UBi 
EIN, | + Pio 3 , rel 5 (Ars) 
> Hpi Pes, Bi | Bi 
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The expected queue lengths are 








Mey 
Fi Ba Ups 
and 
MA... 
Bil 
Eles.) = BIN..] = 

B B 

4 : Bi 


The expected delay times are 


Pon. I 
Fi 
Fi MA ns 
and 
Bi * “MAA : 


The expected waiting times are 


nf 
E(W..] = E[D..J] + — 
loan Fi Uns 
and 
Mei 
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(4.38) 


(4.39) 


(4.40) 


(4.41) 


(4.42) 


(4.43) 





V. COMPARISON OF MODELS 


Comparisons are odorous. 

-William Shakespeare, Much Ado About Nothing 

FORTRAN programs were written to obtain solutions Lou 
the bivariate and decomposition (repairman) models. Results 
for the decomposition (M/M/S) model were obtained from 
previous M/M/S queue calculations made at RAND. Data used 
in the programs was provided by RAND. It consisted of an 
arrival rate, a service rate, and the number of servers at 
each shop of a 20 shop repair facility. It was decided to 
work with only 10 shops of the 20. This decision eliminated 
unnecessary repetitive calculations which would have added 
little to the model comparisons, and still provided enough 
data so that the manpower savings resulting from a consolida- 
tion of back-shop functions could be examined. 

The data could not be directly applied to the models, 
Since arrival rates and service rates were computed from all 
repair transactions at each shop, without differentiation as 
to category of repair, flight-line or back-shop. A rough 
estimate of the ratio of back-shop to flight-line repairs 
(.25) provided a basis for converting the arrival rate into 
the shop into flight-line and back-shop arrival rates. Back- 
shop and flight-line service rates were estimated to be equal, 
SO no conversion was necessary. The number of aircraft used 
in the decomposition models was 25. Table II shows how the 


data was converted. 
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TABLE II 


Data Conversion Table 


Bivariate 







and M/M/S 


Repairman 


Arrival Rate 


Service Rate 


The division of repair personnel between flight-line and 
back-shop portions is not necessary in the bivariate model. 
But the decomposition models require that some division be 
made. For a shop with Ss; personnel, there are Ss; -1 possible 
divisions that may be made. It was decided to make the 
division so that the expected number of customers in the shop 
was minimized. To find this division, the repairman FORTRAN 
Peegtram was run for each shop, with 1 through S; 71 repairmen 


working in the back-shop, and then run again with 1 through 
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S,-1 repairmen working at the flight-line. The expected 


number of customers for each combination of §S and Spar such 


Fi 


. + S.. = S. , was looked at. That combination 
Fi Bl 1 


producing the minimum expected number of customers determined 


that S 


the division of personnel. This division of personnel was 
also used in the M/M/S decomposition model. 

The results of the calculations for the bivariate model 
and the decomposition models are presented in Tables III 
through XII and Figures 8 through 29. In some cases, 
calculations for the back-shop portions of the decomposition 
(M/M/S) model were not made. However, there is still 
adequate information presented to allow a comparison of the 
models. The bivariate model predicts the best performance, 
while the decomposition (repairman) model predicts the 
second-best performance, and the decomposition (M/M/S) model 
predicts the worst performance of the three. The performance 
predicted by decomposition (repairman) model is in most 
cases very close to the performance precicted by the bivariate 
model, with the exception of queue lengths and delay times. 
The bivariate model shows the flexibility of the queueing 
systems presently employed by the USAF, where shop personnel 
can be assigned to the flight-line or the back-shop as they 
are needed. Queueing will not occur until all S; personnel 
are busy in the bivariate model, whereas in a decomposition 


model, queues form when either S ONG Spe personnel are busy 


Sl 
at the flight-line or back-shop of the system under considera- 


tion. The M/M/S queue bad performance is possibly caused by 
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the fact that no reduction in demand for repairmen occurs 


as the number in the system increases. 
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TABLE III 


Measures of Congestion for Shop l 


Calculated from Three Different Models 


Caled) ArLons. Lor 
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= 0.001976, 0.448 
=n OF 010 ac amr = 0.448 
Bivariate Decomposition 
Model Models, Spy = 2,5 Spy = 1 
Measure S; = 3 Repairman M/M/S 
E(Np1) 0.4337 0.4527 0.4647 
Var(Npl) 0.4299 0.4818 N.C. 
E(Qry) 0.001315 0.018/0 0.0227 
‘Var(Qr) 0.00176 0.02702 N.Co 
E(Wpy) Us Laven 203203 2467 
E( Dp) 0,00679 en 9 O72 0.1146 
E(Np1) 0.1098 CieZ 20 Ceo 
Var(Np1) Oly Oss. 0.1361 N.C. 
E(Qp1) 0,0041 0.0129 0.01367 
Var(Qpy) 0.0020 0.0155 N.C, 
E(Wp1 ) 2 BINA 2.4935 2.508 
E(Dp1) 0.0848 0.2614 0.2766 
NeC. - Not calculated 
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TABLE IV 


Measures of Congestion for Shop 2 


Calculated from Three Different Models 


Measure 


E(Np>) 
Var(Np>) 
E(Qp>) 
Var(Qr?) 
E(Wp5) 
E(DF2) 
E(Np2) 
Var (Nx) 
E(Qp9) 
Var (Qp>) 
E(Wp >) 


E(Dp2) 







Calculations» for 
= On0O1 39. UB? =Ooe 2 
= 020056 r Up? = Cay 












Bivariate Deocomposition 








Model 






Models, iy) = 2, SB? = ] 





0.5005 0.5292 0.5467 
0.4970 0, 5760 Meee 
0.00237 0.0290 0.03569 
0.00310 000442 Mer 
3.6938 3.8894 cracks! 
I, Ove 0.2129 0.2563 
0.1280 0.1488 0.1467 
0.1286 0.1636 Nees 
OCs 0.0176 0.0158 
0.0118 0.0219 Nee: 
35 US: 4.184 4.216 
0.339 0.5081 0.5394 


N.C. = Not calculated 
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TABLE V 


Measures of Congestion for Shop 3 


Calculated from Three Different Models 


Measure 


B(NE 3? 
Var (Np3) 
E(QFr3) 
Var(Qr3) 
E(W3) 
E(Dp3? 
E(Ny3) 
Var(Np3) 
E(Qp3) 
Var(Qp3) 
E(W 3) 
E(Dp3) 





Nici. 


33 = 
NS = O2V492, 


Model 


2 hoa 
244804 
0.000006 


0.000008 
UO) e185, 


0.0000059 


wow 97 
0.6/97 
0.000055 
0.000067 
Ug SyaN, 
0.00021 


Bivariate 





oaee 1 t 





0.0123, 
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Calculations for 
Up3 = 0.379 
Up3 = 0.379 












Decomposition 





228811 
US SAE 
0.0097 
0.0176 
2.6475 
0.0089 
0.7886 


0.7691 


0.00154 
0.00211 


2.004 
0.0052 


- Not calculated 


Models, S53 = ve 533 = 4 


Sa 
NoCeo 
0.0469 
N.C. 
206/6/ 
0.0381 
0.8152 
N.C. 
0.0026 
N.C. 
2664/7 
0.0084 
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Beets VL 


Measures of Congestion for Shop 4 


Calculated from Three Different Models 


Calculations 


= 0,00732, Up, = 


BG 
Ary =0,02936 9 Ung = 








pall 


Bivariate Decomposition 
Model Models, Spa = 6 Spq4 = 3 
a> 9 Repairman M/M/S 
E(Nyy) 2.0991 Falaos Ds BGS 
Var(Npy) | 1.9320 1.9962 NeiGe 
E(Qr4) 0.000023 0.0072 Oe02523 
Var(Qry4) | 0.000031 0.0123 N.C. 
E(Wpy,) 3.1949 35 DOE 3 LI 
E(Dpy) 0.000035 0.0107 0.0316 
E(Np,) 0.5235 Oo 52 Nie. 
Var(N,,) 0.5126 me ot N.C. 
E(Qp4) 0.0001 28 0.0040 Nae. 
Var(Qpy) | 0.000179 OS0055 Ne Cs 
E(Wpy) 3.1958 3.2174 N.C. 
E(Dp,) 0.00078 | 0.0225 NG 
Nae. * Not Calculated 
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TABLE VII 


Measures of Congestion for Shop 5 


Calculated from Three Different Models 


Measure 
E(Nps) 
Var(Npo) 
E(Qps) 
Var (Qps) 
E(Wrs) 
E(Dps) 
E(Nps) 
Var(Nxc) 
E(Qgs) 
Var(Q,5) 
E(Wps) 
E(Dg5) 










N.C. 


= Ole 24 9 Ups = 


Calculations for 








= 0,0438 9 Urs 


Bivariate 


Model 
ss. = eZ 


Zia Zs 
PaDUZ9 
0,.00000084 


0,00000107 
2.6810 
0.00000079 
0./0245 
0.682/ 
0.00269 
0.0524 
266811 
0.01039 


Repairman 
Peo 943 
DS TAGS) 
0.00214 
0.0035 
226830 
0.0020 
Oe /e59 
OPT es 
0.0016 
0.0022 


2.63864 
0.00541 


- Not calculated 
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Os Sis 


Decomposition 


Models, Sps = 8, Sps5 = 4 


M/M/S 
Saftey 
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228166 
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TABLE VIII 


Measures of Congestion for Shop 6 


Calculated from Three Different Models 







Calculations for 
SOG? Aly, Up6 = On 7S 
Oe 75 

















Bivariate Decomposition 









Model 






Models, SEr6 = var SS tm 3 





Measure | 


S$. = 10 M/M/S 
E(Np6) 2.545 2.616 3,000 
Var(Npe) | 2.286 2.3661 N.Co 
E@e) 0.00001843 | 0.0051 0.0839 
Var(Qpe) | 0.000025 0.00875 Ngee 
E(Wp¢) 1.290 1.2928 1.3274 
E(Dre) 0.00000935 | 0.0025 0.0371 
E(Npe) 0.635 0.7161 N.C. 
Var(N,-) | 0.619 0.7207 NeCo 
E(Qxe) 0.000138 0.00922 NAO 
Var(Qpe) | 0.000182 0.01365 N.C. 
E(Wpe) 1.2903 1.3072 N.C. 
E(Dpe) 0.00028 0.0168 RECs 





NeC.e = Not calculated 


77 





PROBABILITY OF N AIRCRAFT DOWN 


oO 


29 


20 


PLOT OF PROBABILITY VS. AIRCRAFT DOWN 
FOR FLIGHT-LINE REPAIR AT SHOP NR 6 


LEGEND 

O BIVARIATE MODEL 

A REPAIRMAN MODEL 
Rg ry M/M/S MODEL 


>O 


OD 


Gl 
wo 
& 2 
oO l 2 3 A 5 S 7 gS 
N 
FIGURE 18 


78 





PROBABILITY OF N AIRCRAFT DOWN 





PLOT OF PROBABILITY VS. AIRCRAFT DOWN 
FOR BACK- SHOP REPAIR AT SHOP NR 6 


LEGEND 
O BIVARIATE MODEL 
A REPAIRMAN MODEL 


FIGURE 19 


ges 





TABLED x 


Measures of Congestion for Shop / 


Calculated from Three Different Models 


Calculations for 


B/ 


r7 = 000368 ,'p7 = 0.298 





Bivariate Decomposition 


Model Models, Sp7z7 = 4y Sp7 

E(Npo) 246894 341859 4» 936 

Var (Np7) 204881 Ba \) SS N.C. 
E(Qp7) 0,0257 0.49206 1.8490 
Var(Qp7) | 0.0493 1.3335 NeC. 
E(Wp7) 3.3880 Bs SSS Be o054 
E(Dr7) 0.0323 0.6130 2.0097 
E(Np7) 0.7401 0.8476 0.9069 
Var(Np7z) | 0.7725 1.0286 N.C. 
E(Q37) 0.0756 0.1019 0.1351 
Var (Qp7) G21 335 0.1869 iN Cs 
E(Wp7) 3.7297 3.8145 3.943 

E(Dp7) 0.3310 0.4587 0.5872 





NoCe = Not calculated 
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TABLE X 


Measures of Congestion for Shop 8 


Calculated from Three Different Models 





Calculations for 








Ce00s 2, Hog = Ob4y 2 






0.4/2 





Dat 
bry 
ee) 

| 





Bivariate Decomposition 















Model Models, Spe = =) Spg = 2 




























E(Npg) 151426 1.1570 Ih Zale 
Var (Nrg) 1.0904 Lugs 7 NoCo 
E(Qrg) 0.000013 0.00126 0.0027 
Var(Qrg) | 9.000017 0.0018 NeCo 
E(Wre) DNS 2.1210 Dio 4 
E(Dzg) 0.000025 0.0023 0.0048 
E(Ngg) 0, 2857 0.3053 Nees 
Var(Npe) G5 ES 2S 0.3135 Paes 
E(Qpg) 0.000051 0.0060 NC. 
Var(Qpeg) | 0.000069 0.0078 N.Co 
E (Wag) Ze Mois Pesnoiles N.C. 





E(Dgg) 0.000378 0.0426 NaS 


NoCe = Not calculated 


8 3 





PROBABILITY OF N AIRCRAFT DOWN 


40 PLOT OF PROBABILITY VS. AIRCRAFT DOWN 


FOR FLIGHT-LINE REPAIR AT SHOP NR & 


® 
0) 
i) 
LEGEND 
© BIVARIATE MODEL 
A REPAIRMAN MODEL 
30 & cl M/m/S MODEL 
.Q0 
a 
8 
20 
1S 
PO 
a 
3) 
05 
B 
© 
0 < () Cs) {-) t:) 
0 } 2 3 4 5 6 7 & 
’ N 
FIGURE 22 


84 





PROBABILITY OF N AIRCRAFT OOWN 





PLOT OF PROBABILITY VS. AIRCRAFT DOWN 
FOR BACK-SHOP REPAIR AT SHOP NR 8B 


LEGEND 
© BIVARIATE MODEL 
A REPAIRMAN MODEL 
cm M/M/S MODEL 


FIGURE 23 


85 





TABLE XI 


Measures of Congestion for Shop 9 


Calculated from Three Different Models 


Measure 


E(NF9) 
Var (N29) 
E(Qr9) 
Var(Qrq) 
E(Wro) 
E(Drg) 


E(Np 9) 
Var(Npo) 
E(Qgo) 
Var(Qp9) 
E(Wp9) 
E(Dg9) 





Model 


4.93824 
309394 
0.000000 
0.000000 
3,003 
0,000000 
1.2490 
1.1866 
0.000000 
0.000000 
35 OST 
0.000000 


Calculations 
An9 = One he. 
A. = 0,0884 r 


Bivariate 





Sg = ee 










D0 2444 
4.1442 
0.000000 
0.00001 
3.003 
0.000000 
1.5599 
1.4626 
6 HON? 
0.00002 
CPuelek: 
0.00003 





NeCe- Not calculated 


§6 


UBS 
Ung 


Decomposition 





— O83 3 


Oyeo0 





Models, Seg = 15, S39 a 


Repairman 


M/M/S 


6.639/ 
NoCo 
0.0031 
NoCo 
3.004 
0.0014 
NoCo 
N.C. 
NeCe 
N.Co 
N.C. 


NeCo 





PROBABILITY OF N AIRCRAFT DOWN 


.30 


PLOT OF PROBABILITY VS. AIRCRAFT DOWN 
FoR FLIGHT-LINE REPAIR AT SHOP NR Q 


LEGEND 

@ BIVARIATE MODEL 
A REPAIRMAN MODEL 
O M/M/S MODEL 


FIGURE 24 


87 








PROBABILITY OF N AIRCRAFT DOWN 


,35 PLOT OF PROBABILITY vS. AIRCRAET DOWN 
FOR BACK-SHOP REPAIR AT SHOP NR 9 


4 

30 LEGEND 

O BIVARIATE MODEL 

@ A REPAIRMAN MODEL 
A 
(en 
Oo 
20 & 
15 
rN 
10 
© 
05 4 
© 

yen 

© 
O & 

eo) 2 3 4 5 6 i ee 
N 
FIGURE 25 


88 





TABLE XII 


Measures of Congestion for Shop 10 


Calculated from Three Different Models 


Galculatlonseror 


10 0.00532 ¥ a) (he 0,389 


Bivariate Decomposition 


Model Models, Spy] 0=2 ’ Sp o=l 
Measure 





S19 = 3 M/M/S 
E(NF10) 0.3366 0.3458 tes 522 
Var(NF10) | 0.3335 Wo Sz Noa 
E(QF10) 0.00005 0.00865 0.01029 
Var(QFr1 9) | 0.00006 CeO 155 NoCo 
E(Wr1 0) De SHES 2.6366 2.648 
E(Dp10) 0.0038 0.0659 0.0774 
E(Npi0) 0.0346 0.0926 0.0896 
Var(NpjQ) | 0.0845 0.1003 Nin Oc 
E(Qz10) 0.000712 0.00753 0.00737 
Var(QpjQ) | 02000307 0.00874 NoCo 
E(Wp10) 2.592 207982 2.801 





E(Dgi9) | 020218 0.22759 0.2304 


NeCo - Not calculated 
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VI. CALCULATION OF MANPOWER SAVINGS 


When wilt thou save the people? 
Oh, God of Mercy! when? 
The people, Lord, the people! 
Not thrones and crowns, but men! 
pueenezes EH11e6t, Poetical Works 

One application of the models developed in the previous 
chapters is the determination of manpower savings attainable 
from a scaled-up facility. Specifically, the savings 
achievable if the back-shop functions of five identical bases 
were consolidated into a large back-shop facility was deter- 
Mined. Each base was assumed to consist of 10 shops, each 
identical to Shops 1 through 10 used in the model comparisons 
made in Chapter V. Each base operates 25 aircraft with 
failure characteristics to those used in the model comparisons. 
The scheme was to leave crews at the 5 bases to handle only 
flight-line repairs. Back-shop repairs were made at the 
consolidated back-shop. 

The decomposition (repairman) model was used to determine 
the number of repair personnel to leave at each base for 
flight-line duties. The back-shops were manned from the pooi 
of remaining maintenance personnel. The consolidated back- 
shops were modeled as M/M/S queues since the number of 
aircraft failing was large enough (125) so that the model 
predicts reasonably well. 

Each base was manned with a total of 87 repair personnel 


at the ten shops, with 58 assigned to flight-line duties and 
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29 assigned to back-shop duties (according to the repairman 
model). For 5 bases there were 435 total men available for 
the base consolidation, with 290 assigned to flight-line 
duties, and a pool of 145 men available for assignment to 
the consolidated back-shops. 

It was decided to assign men to the consolidated back- 
shops so that the expected number of aircraft down at each 
shop was within 10% of the number predicted by the decomposi- 
tion (repairman) model, with shops manned as in Chapter V. 
For example, if the expected number of aircraft in the back- 
shop of Shop i was 1.00 before the consolidation, then 
5 bases would have, on the average, 5.00 men in the back-shops 
of Shop i. The standard for the consolidated back-shop would 
then be an expected number of aircraft down of less than 5.5. 

The arrival rate used at the consolidated back-shop of 
Shop 1 was 125A,,- The service rates were unchanged. The 
M/M/S queue model was run for various manning levels, until 
the expected number of aircraft in the shop dropped below 
the standard discussed above. Table XIII presents the result 
of the calculations of four measures of congestion made for 
the consolidated back-shop facility, and also shows the 
manpower savings achieved. 

A total of 89 personnel were eliminated by consolidation 
of back-shop repairs. This is roughly the former manning 
level of one base, or a savings of roughly 20% achieved by 
scaling-up the back-shop functions. Recall, that the total 


number of shops at the base from which the data was collected 
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was twenty. The number of men saved by reorganization should 
be roughly twice that predicted here, or a little less than 


200 repair personnel. 
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VII. CONCLUSIONS 


‘Is there any point to which you wish to draw 
my attention?' 

‘To the curious incident of the dog in the night- 
time.' 

"The dog did nothing in the night-time. ' 

‘That was the curious incident,' remarked Sherlock 
Holmes. 


The decomposition (repairman) model was in close 
agreement with the bivariate model, except with respect 
to queue lengths and delay times. It is quite suitable for 
use in problems such as base consolidation studies, where 
Savings of repair personnel are sought. The decomposition 
(M/M/S) model provides conservative predictions, which makes 
it more suitable for studies made for a risk-averse decision 
maker. 

If performance predictions are desired for repair 
configurations similar to that described in Chapter II, 
where the effects of sustained changes in tempo of operations 
or the effect on system performance of manpower reduction or 
augmentation is desired, then the bivariate model is more 
appropriate than either of the two decomposition models. 
Furthermore, although all three models are over-simplified, 
the balance equations describing more complex situations can 
easily be written down. Then they may also be solved using 
Gauss-Seidel iteration. Thus techniques to replace complex 
Simulation, at least in some instances, are at hand. Such 


work is under continuing development. 
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